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: ABSTRACT 

^ ■ Results are presented from a programme of detailed longslit spectroscopic observa- 

Psj . tions of the extended emission-line region (EELR) associated with the powerful radio 

^ \ galaxy PKS 2356—61. The observations have been used to construct spectroscopic dat- 

\ acubes, which yield detailed information on the spatial variations of emission-line ratios 

G> ■ across the EELR, together with its kinematic structure. We present an extensive com- 

^ . parison between the data and results obtained from the MAPPINGS II shock ionization 

code, and show that the physical properties of the line-emitting gas, including its ion- 
O \ ization, excitation, dynamics and overall energy budget, are entirely consistent with a 

^ • scenario involving auto-ionizing shocks as the dominant ionization mechanism. This has 

ci ■ the advantage of accounting for the observed EELR properties by means of a single phys- 

j> I ical process, thereby requiring less free parameters than the alternative scheme involving 

photoionization by radiation from the active nucleus. Finally, possible mechanisms of 
shock formation are considered in the context of the dynamics and origin of the gas, 
specifically scenarios involving infall or accretion of gas during an interaction between 
the host radio galaxy and a companion galaxy. 

Subject headings: galaxies: active — galaxies: individual (PKS 2356—61) — galaxies: 
ISM — galaxies: kinematics and dynamics — radio continuum: galaxies — shock waves 
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1. Introduction 

Powerful radio galaxies are often observed to con- 
tain systems of ionized gas, extended up to scales of 
several tens of kpc from the nucleus. Detailed studies 
of these extended emission-line regions (EELRs) have 
revealed a wide range of morphological and kinematic 
properties (e.g., Baum & Hcckman 1989a, 1989b; 
Tadhunter, Fosbury & Quinn 1989; McCarthy, Spin- 
rad, & van Breugel 1995): in some objects the EELR 
is coincident with the radio lobes, while in the major- 
ity of cases the line-emitting gas is extended at modest 
angles away from the radio axis and does not display 
a direct spatial association with the radio emission. 
The EELR kinematics (Tadhunter et al. 1989; Heck- 
man et al. 1989; Baum et al. 1992) can be gener- 
ally classified into "rotators" and "non-rotators" , the 
latter being further subdivided into "calm" and "vi- 
olent". The rotator EELRs typically exhibit circu- 
lar velocity fields approximately consistent with the 
massive halos inferred for radio ellipticals, while the 
non-rotators display less well-ordered kinematics and 
little or no evidence of settled dynamics. 

Explanations for the origin of the gas include merg- 
ers or interactions involving gas-rich galaxies; large 
differences between the gaseous and stellar kinemat- 
ics (Phillips et al. 1986; Tadhunter et al. 1989; Bet- 
toni et al. 1990) suggest external infall of gaseous 
material on recent timcscales compared to the age 
of the elliptical (Habe & Ikeuchi 1985, 1988; Varnas 
et al. 1987). Detailed kinematic studies of EELRs 
often reveal a correspondence between the radio axis 
and the gas "rotation" axes (e.g., Simkin 1977; Heck- 
man et al. 1985; Tadhunter et al. 1989), suggesting 
a direct connection between gas infall and nuclear ac- 
tivity. It has also been suggested that the gas has 
condensed out of a hot intergalactic medium (Fabian 
et al. 1981; Canizarcs et al. 1987); however, the an- 
gular momentum of gas in radio galaxies tends to be 
higher than that of typical cooling flow EELRs. Fur- 
thermore, the importance of mergers in fuelling ga- 
laxy activity is suggested by the typical morphologi- 
cal disruption of host galaxies (Hutchings 1987; Smith 
& Heckman 1989) as well as apparent environmental 
differences between radio-loud and normal ellipticals 
(Lilly & Longair 1984; Spinrad & Djorgovski 1984). 

Thus, a principal question in the study of the for- 
mation and evolution of radio galaxies concerns the 
detailed physical properties of the EELRs. In particu- 
lar, models of the EELR ionization and dynamics are 



of importance in understanding their origin, as well 
as their role as a source of fuel for the active nucleus, 
related to the more general question of the triggering 
mechanism for active galaxies. 

Excitation of radio galaxy EELRs has generally 
been considered in the context of "unified schemes" 
of radio galaxies and quasars (e.g., Barthel 1989), in- 
volving photoionization of the gas by radiation from 
the active nucleus. This is analogous to the sce- 
nario inferred on kpc scales for EELRs in Seyfert 
galaxies, where nuclear photoionization is supported 
by "ionization cones" in line-ratio maps in a num- 
ber of objects, and by the presence of polarized 
broad emission lines (e.g., Antonucci & Miller 1985; 
Mulchaey, Wilson, & Tsvetanov 1996). While evi- 
dence for orientation-based unification of radio ga- 
laxies and quasars has been inferred from a range 
of observational studies, including radio galaxy mor- 
phology and beaming studies. X-ray properties, and 
host galaxy populations (e.g., Antonucci 1993 and ref- 
erences therein; see also Gopal-Krishna, Kulkarni, & 
Wiita 1996), it is less clear that all the detailed prop- 
erties of the line-emitting gas in all radio galaxies 
can be explained in terms of nuclear photoionization. 
Firstly, radio galaxy EELRs are extended on much 
larger scales than those in Seyfert galaxies, typically 
by at least an order of magnitude, and do not display 
clear ionization c;onc;s. Furthermore, relationships be- 
tween EELR luminosity and core power or total ra- 
dio power typically display a scatter of 1 — 2 orders 
of magnitude (e.g., Baum & Heckman 1989b), which 
may be due in part to different ionization mechanisms 
dominating in different objects that have otherwise 
similar core properties. Another interesting effect is 
that the global kinematic properties of EELRs have 
been shown to correlate with their overall excitation 
state; this is not a directly expected consequence of 
unified schemes, and suggests the importance of lo- 
cal ionization in some classes of objects (Baum et al. 
1990, 1992). 

Moreover, detailed studies of individual radio ga- 
laxies have shown that the EELR properties cannot 
always be readily accounted for in terms of nuclear 
photoionization. Specifically, the extended material 
in a rnimbcr of objects display emission-line fluxes 
and excitation states that are inconsistent with the 
nuclear ionization properties required to photoion- 
ize the material closer to the AGN (e.g., Danziger 
et al. 1984; van Breugel et al. 1986; Tadhunter et al. 
1988; Meisenheimer & Hippelein 1992). In such ob- 
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jects, the importance of local ionizing mechanisms, 
such as shocks, has been postulated in order to ex- 
plain the large-scale properties of the extended gas. 
Thus, while orientation effects are clearly likely to af- 
fect some observational properties of radio galaxies 
and quasars, we point out that the ionization and ex- 
citation properties of the EELRs may be dominated 
to varying degrees in individual objects by local ion- 
ization mechanisms, related for example to shocks. 

In this paper we address in detail the feasibility of 
"auto-ionizing shocks" as a viable ionization mecha- 
nism for the extended gas. One possible origin for 
shocks is collisions between streams of gas during a 
galaxy merger. This scenario is motivated by the 
"tidal" EELR morphology often observed in radio 
sources, where the gas is neither conical nor directly 
associated with the radio plasma. In objects where 
the radio plasma coincides with the EELR (e.g., Pri- 
eto et al. 1993; Tadhunter et al. 1994) it is possible 
that associated bowshocks may play a direct role in 
ionizing the gas. For example, the auto-ionizing shock 
model has been applied successfully to observations of 
emission-line gas in Centaurus A (Sutherland, Bick- 
nell & Dopita 1993). The physical appeal of auto- 
ionizing shocks is that the gas energetics, excitation 
and kinematics are accounted for in a single, physi- 
cally self-consistent model, requiring less free param- 
eters than photoionization by a hidden AGN. 

2. Observational Programme 
2.1. Target Object 

In order to obtain as complete a description as pos- 
sible of the physical properties of the gas, a combi- 
nation of kinematic and excitation information is re- 
quired across its entire extent. We describe the results 
from a program of low- and high-dispersion longslit 
spectroscopy of the radio galaxy PKS 2356—61, sam- 
pling the EELR completely through the use of multi- 
ple slit positions. 

The radio galaxy PKS 2356—61 is one of the 
strongest southern FR II sources, with a total radio 
power Pi.4GHz ~ lO^''^'^ WHz""'^, and was discovered 
during the "Mills Cross" radio survey of the southern 
sky by Mills et al. (1961). Its double-lobed mor- 
phology has been apparent since the first radio syn- 
thesis observations were obtained (Ekers 1969). The 
optical counterpart was identified by Westerlund & 
Smith (1966) to be a type E3 or D elliptical of V 
magnitude 16, in a sparsely populated field (contain- 



ing 10 other galaxies within a region ~ 21' in size, i.e., 
~ 2.3 Mpc for {Ho = 75kms-iMpc-\ qo = 0, which 
we use throughout this paper). The first redshift de- 
termination was by Whiteoak (1972), who measured a 
value of z = 0.0959 based on optical emission-lines. A 
nuclear spectrum obtained by Danziger & Goss (1983) 
shows strong, high-excitation narrow-line emission, 
with a high [OIII]A5007/H/3 ratio (^ 10) and the pres- 
ence of significant [NeV]A3426. This is also evident 
in the fluxes published by Robinson et al. (1987) for 
the nuclear region and two off-nuclear regions (7" E 
and 4" W). These and other optical properties are 
summarized in Tabic |^. 

The only previously published emission-line image 
is in Ha (Danziger & Focardi 1988), showing a bright, 
elongated emission-line region centred on the galaxy 
nucleus, with an arc or shell- like feature on its eastern 
side. A very faint tail extends from the east to the 
south; a high-dispersion longslit spectrum of one of 
the blobs in this tail (Tadhunter et al. 1989) shows a 
relatively small velocity gradient (< 200kms~^). 

2.2. Observations and Data Reduction 

The observations were obtained during dark time 
on the nights of 1 and 2 August, 1992, using the RGO 
Spectrograph mounted at the f/8 Cassegrain focus of 
the Anglo- Australian Telescope ( AAT) . The detector 
used was the Tektronix #2 CCD, which consists of 
1024 X 1024 24^m pixels, each corresponding to 0'.'77 
on the sky. Details of the observations are presented 
in Table ^. Conditions were sufficiently photometric 
that reliable absolute photometry was obtained, and 
the seeing ranged between 1'.'2 and 1'.'5. A detailed 
account of the observations and reduction procedures 
is presented in Koekemoer (1996). 

Longslit exposures were obtained at parallel posi- 
tions across the entire extent of the EELR, thereby 
producing two spectral datacubes (at high- and low- 
dispersion) in Right Ascension, Declination and wave- 
length. New slit positions were acquired by offsetting 
the telescope in increments of 1". This is most reliably 
achieved through the use of "Star-to-Probe" guiding: 
the guiding assembly is moved by a specified amount, 
and the telescope moves automatically in order to re- 
position the guide star on a standard location. By 
re-acquisition of the galaxy nucleus and nearby stars 
throughout the night, it was verified that positioning 
errors were <^ 0'.'15. 

Longslit observations are susceptible to problems 
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Table 1 

General properties of PKS 2356-61. 



Notes 



R.A. (B1950) 


23'^ 56" 29!37 ± 0!02 


1,2 


Dec. (B1950) 


-61° 11' 40'.'6 ± C/.'S 


1,2 




0.0963 


3 




0.0959 ±0.0001 


3,4 


'^abs+em 


0.0958 ±0.0003 


5 


Distance Dl 2:abs+em (Mpc)^ 


383 




Scale (kpc/arcsec)'f 


1.86 




Galaxy type 


D/E3 


1,2,6 


V mag. 


16 


1 


B mag. 


17 


6,7 


Continuum extent 


6" X 4" 


1 


Continuum P. A. 


6° 


1 


EELR extent 


14" X 8" 


8 


-P'[OIII]A5007/ -fH/3 

FH/3(erg/s/cm2) 
A\4ot (kms-i) 


13.2 


9 


3.94 X 10-15 


9 


180 


8 


AKiax (kms^^) 


90 


8 


Pl.4GHz (WHz-1) 


6.3 X 10^5 


10 



Note.— 1: Westerlund & Smith 1966; 2: Ekers 1970; 3: Tritton 
1972; 4: Whiteoak 1972; 5: Danziger & Goss 1983; 6: Sutton 1968; 
7: Tadhuntcr ct al. 1989; 8: Danziger & Focardi 1988; 9: Robinson 
et al. 1987; 10: Wright & Otrupcek 1990. 
+ We use Ho = 75kms-i Mpc-\ qo = 0. 
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Table 2 

Spectrophotometric Observations of PKS 2356-61. 





1 August 1992 


2 August 1992 


Grating 


270R 


1200V 


Wavelength/pixel 


3.39 A 


0.58A 


Spectral range 


3996-^74651 


5305 ^5895 A 


Emission lines 


[OII]A3727^[SII]A6731 


H/3, [OIII]A5007 


Calibration arcs 


He/Ne/Ar 


Cu/Ar 



Note. — The data obtained each night consisted of exposures 
in 14 parallel, adjacent slit positions at a position angle of 0°. 
The slit width was 1" and the slit positions were offset in R.A. in 
increments of 1"; the spatial coverage extended from 5" W ^ 8" E 
of the nucleus. The exposure time at each slit position totaled 
1200 sec. 



due to atmospheric effects. Seeing fluctuations can 
cause changes in the contribution from regions imme- 
diately outside the slit, while low-dispersion data can 
suffer from differential refraction, resulting in vary- 
ing contributions from regions outside the slit as a 
function of wavelength. However, all the data were 
obtained at zenith distances ^ 30 — 40°, thus dif- 
ferential refraction effects are <^ 0'.'2 along the wave- 
length regime covered. The impact of these effects 
was further minimized by employing a slit width of 
1", comparable to the seeing disk size, and by spa- 
tially smoothing the final dataset using 3 x 3-pixel av- 
eraging. 

Smooth-spectrum stars and photometric standards 
(Baldwin & Stone 1984) were observed to correct for 
atmospheric absorption lines and the detector/gra- 
ting spectral response. Routines in the IRAF "one- 
dspec" package were used to apply these corrections 
to the data. Frames were de-biased using the cor- 
responding overscan regions, and no two-dimensional 
bias structure was present on the detector. Pixel-to- 
pixel flat-fields were obtained from exposures of the 
internal tungsten lamp. The frames were corrected for 
optical vignetting along the slit using twilight expo- 
sures obtained with the same slit width. The spectra 
were corrected for atmospheric extinction using stan- 
dard extinction curves for Siding Spring Observatory 
(available within IRAF). 



Internal arc exposures were obtained at regular in- 
tervals throughout each night, in order to correct for 
instrumental shifts. A two-dimensional wavelength 
transformation for each arc frame, subsequently ap- 
plied to the data, was obtained by fitting arc lines at 
multiple positions along the slit, using routines in the 
IRAF "twodspec" package. The resulting residuals in 
wavelength were less than 0.1 pixels across the regions 
of interest. After wavelength calibration, the night- 
sky spectrum for each frame, obtained from empty 
regions along the slit, was subtracted from the data. 

3. Results 

3.1. Morphology 

We used the low-dispersion datacube to investigate 
the spatial distribution of the stellar continuum and 
emission-line gas in PKS 2356—61. The datacube 
was smoothed in R.A. and Dec. using a "boxcar" 
filter (3x3 pixels), in order to minimise effects due 
to possible changes in seeing throughout the night. 
A continuum image, shown in Figure |l|a, was pro- 
duced by summing all channels with observed wave- 
lengths in the range 5591 - 6685A ( i.e., rest wave- 
lengths of 5100— 6IOOA, corresponding approximately 
to y-band), taking care to avoid the 5577A night sky 
line as well as [OIIIJA5007 and other strong emission 
lines. We fitted elliptical isophotes to the host ga- 
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laxy in this image using the IRAF/STSDAS package 
"elUpse" (Jedrzejewski 1987), and compared the re- 
sults with ellipses fitted on a broad-band CCD image 
of the object. In neither case did we find any sig- 
nificant residual continuum features co-incident with 
the emission-line region; instead, the stellar contin- 
uum follows a well-behaved elliptical light distribu- 
tion. Furthermore, the radial distribution of sur- 
face brightness, ellipticity and position angle as de- 
rived from the spectroscopic datacube showed excel- 
lent agreement with those obtained from the CCD im- 
age, confirming that the datacube is capable of yield- 
ing spectrophotometric results as reliable as those ob- 
tained using conventional imaging photometric tech- 
niques. 

Some faint, extended continuum emission appears 
to be associated with the objects toward the north 
of the source; this is also evident on the broad-band 
CCD image. Spectroscopically, this structure is as- 
sociated with Ha emission at very similar systemic 
velocities to the gas in the central source; it is there- 
fore likely that this material is dynamically associated 
with the central galaxy. The brighter objects, how- 
ever, are foreground stars, as confirmed by examina- 
tion of their summed low-dispersion spectra. 

"Narrow-band" images of the emission-line gas in 
[OIIIJA5007 and Ha are presented in Figures |l|6,c. 
These were obtained by parameterizing the line pro- 
files in terms of Gaussian components. All pixels con- 
taining detectable line emission were fitted using the 
"specfit" routine within IRAF, yielding values and er- 
ror estimates for the fluxes, velocities and widths of 
the lines [OHI]A5007, H/3, Ha, [NH]AA6548,6583 and 
[SH]AA6716,6731. The spectral resolution is - 4A, 
sufficiently high to separate the Ha and [NH] pro- 
files but also low enough that each line profile can be 
adequately described by a single Gaussian. The re- 
sulting images in Figures |l|6,c represent all pixels with 
flux values greater than the corresponding formal 3-a 
limits obtained from the fits. 

The emission-line region associated with the host 
galaxy may be described generally in terms of two 
components, a bright elongated central region and a 
fainter tail extending up to 15 — 20" from the nu- 
cleus. The tail itself consists of three distinct clouds 
which we designate A, B and C, in order of increasing 
distance from the central region. No association with 
the radio axis (at PA ~ 135°) is evident, either in the 
central region or in the tail. However, it is interest- 
ing to note that the major axis of the central region 



is approximately perpendicular to the orientation of 
the host galaxy optical continuum distribution. 

A large Ha-emitting region is evident to the north 
of the host galaxy, displaying relatively little [OIII] 
A5007 emission. This region corresponds to the faint, 
diffuse emission apparent in the continuum image. A 
faint extension can be seen on cither side of its cen- 
tral component; the total projected length of these 
extensions is 55 kpc. The systemic velocity of this 
object is within ^ 200kms~^ of that of the central 
galaxy, and this is considerably smaller than expected 
if two cluster objects were randomly superimposed (in 
which case the velocity difference should be compa- 
rable to the projected cluster velocity dispersion, i.e., 
^ 1000 km s^^). Thus, it is possible that the Ha ob- 
ject corresponds to a gas-rich galaxy that is interact- 
ing with the host galaxy of PKS 2356—61, in which 
case its extended morphology may be accounted for 
by tidal disruption. 

3.2. Kinematics 

The high-dispersion spectral datacube has a veloc- 
ity FWHM resolution ^ 70kms~^, the spectral sam- 
pling being 0.8A. The observations were taken at the 
same slit positions and offsets as those used for the 
low-dispersion datacube, but the data were obtained 
on a different night, so that seeing effects could po- 
tentially introduce differences in the flux distribution 
compared to the low-dispersion datacube. In order to 
minimize such effects, the dataset discussed here has 
been smoothed in the same way as the low-dispersion 
dataset, by using a "boxcar" filter of 3 x 3 pixels. 

In Figure |^ we present several different views of the 
continuum-subtracted [OIIIJA5007 datacube. These 
images were created with the AVS package]^. The 
orientations in all panels are such that redshift in- 
creases towards the right. The top left panel (Fig. ||a) 
shows the flux distribution of the gas as viewed on 
the sky, with north toward the top and east to the 
left. The cube is then presented at a number of differ- 
ent viewing angles; for example. Figure ^fl shows the 
equivalent of a longslit spectrum covering the whole 
emission-line region — in this image, north is at the 
top and south is at the bottom (i.e., the slit direction 
runs vertically along the page). 

A more quantitative representation of the velocity 
structure is shown in Figure 0; a comparison of these 
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Fig. 1. — Images of PKS 2356—61 derived from the low-dispersion spectral datacube. Positions are given in offsets 
of R.A. and Dec. relative to the nucleus; North is to the top and East is to the left in all images, (a) "F-band" 
continuum image; (b) [O111JA5007 image; (c) Ha image. In all images the slit is oriented north-south and pixel sizes 
are 1" x O'.'TT (corresponding to a slit width of 1", with a spatial scale along the slit of C'TT/pixel). Pixels along the 
slit have not been rebinned into 1" increments, in order to avoid the introduction of artifacts due to non-uniform 
resampling of flux. Instead, the aspect ratio of pixels is plotted in such as way that the spatial scale in arcseconds 
is the same along each axis. 
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Fig. 3. — High-dispersion [OIIIJA5007 velocity structure of the emission-line gas in PKS 2356—61. Each spectrum 
represents the sum of a region ~ 3 — 4" across, i.e., larger than the scale used for smoothing the datacube, thus 
adjacent spectra are independent of one another. All velocities are relative to the systemic velocity (mean emission- 
line redshift of the entire system), which is indicated by the vertical axis. The orientation of this image is the same 
as in Figure H; North is toward the top and East is to the left. 
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data with Figure g reveals that the three clouds com- 
prising the tail are spatially and kinematically dis- 
tinct, having different central velocities as well as dif- 
ferent internal velocity fields. However, the velocity 
widths of individual components in these clouds are 
almost all the same, ranging from 120 to 150kms~^. 
Here we describe the general characteristics of the tail 
clouds. 

Cloud C has a higher redshift than any of the other 
gas components in the EELR. It contains only one 
velocity component, with an almost constant veloc- 
ity width. The velocity centroids also show little or 
no systematic change across the entire extent of this 
cloud. 

Cloud i? is at a lower redshift and is comprised of 
two kinematically distinct components, Bl and B2, at 
lower and higher redshifts respectively. These compo- 
nents both have a very similar velocity width to the 
gas in C, and the width also appears to remain con- 
stant across their entire extent. Furthermore, they 
overlap almost entirely on the sky and their surface 
brightnesses are very similar. Component Bl shows 
no systematic change in velocity centroid but B2 dis- 
plays a smooth velocity gradient across its entire ex- 
tent, being more redshifted in the direction of C. A 
faint connection may be seen between B2 and C, but 
this is possibly a result of the smoothing since the un- 
smoothed dataset contains a set of pixels between the 
two clouds where there is no detectable [OIH]A5007 
emission. The velocity gradient is not a result of 
the smoothing since it extends over many smoothing 
scales. 

Cloud A contains a bright component with a width 
similar to that of cloud Cbut it also contains a fainter, 
broader, blue-shifted component. The bright compo- 
nent shows no systematic change in its velocity cen- 
troid, which is intermediate between those of the two 
components in B. Its velocity centroid is also very 
close to the value of the systemic velocity. The faint 
broad blue extension in A is characteristic of a "blue 
wing" although given its low S/N, it could also be a 
single separate component, as is the case for B. The 
faint component is considerably broader, having ap- 
proximately the same width as the total difference 
between the two components in B. A significant spa- 
tial overlap is apparent between clouds A and B, even 
in the unsmoothed datacube. 

The central component is kinematically distinct 
from the clouds in the tail. It is dominated by a broad 
component, with FWHM values 300 - 400kms~^, 



and also has a much higher surface brightness than 
any of the tail regions. This region also contains 
asymmetric extensions to the main component, pos- 
sibly indicating narrow components similar to those 
found in the tail clouds; these could also be due to 
the presence of blue and red wings. 

3.3. Emission-Line Ratios 

We next consider the physical conditions and exci- 
tation state of the gas, using the spatial variation of a 
number of emission- line ratios across the EELR as de- 
rived from the low-dispersion datacube. In this sub- 
section we present several emission-line ratio maps, 
followed by composite spectra for a number of sub- 
regions within the EELR. 

3.3.1. Line- Ratio Images 

In Figure ^ we present images of the spatial varia- 
tion of the line ratios [OIII]A5007/H/3, [NII]A6583/ 
Ha, [SII]AA6716-f6731/Ha, [SII]AA6716/6731, and 
Ha/H^. 

An important consideration in quantifying the Bal- 
mer line fluxes is the degree of Balmer absorption 
present in the underlying stellar continuum. This 
has been addressed in the literature using a variety of 
techniques (e.g., Koski & Osterbrock 1976; Filippenko 
& Sargent 1988; Storchi-Bergmann et al. 1992), gen- 
erally involving an independent determination of the 
stellar continuum. However, in PKS 2356—61 there 
are no strong emission-line free regions of continuum, 
and the continuum is also too faint to allow an unam- 
biguous determination of the intrinsic degree of ab- 
sorption, for example by using synthetic stellar spec- 
tra. 

Instead, the significance of stellar absorption may 
be determined through a comparison of the variation 
of the ratio Hq;/H/3 and the stellar continmmr. It can 
be seen in Figure ^ that a strong peak in this ra- 
tio is present at a location immediately to the east 
of the nucleus, increasing by almost a factor of 2 to 
values in the range 5 — 6. Since the stellar continuum 
is relatively weak, a drastic change in stellar popula- 
tion would be required to account for this effect by 
Balmer line absorption; however, no such change is 
observed in the host galaxy continuum spectra. Thus 
we infer that the variations in Ha /H/3 across the cen- 
tral EELR are indicative of intrinsic variations in the 
amount of extinction, and not the result of underlying 
absorption. In the tail, however, there is a possibil- 
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ity of significant underlying absorption in the contin- 
uum, which we discuss further when presenting the 
one-dimensional low-dispersion summed spectra. 

Examination of the line-ratio images reveals the 
presence of a "band" of excitation across the nu- 
cleus, extending towards the south-west and north- 
east. This feature is particularly apparent on the im- 
ages of [OIII]A5007/H/? and [NII]A6583/Ha, display- 
ing typical line ratio values of [OIII]A5007/H/? 10, 
[NII]A6583/Ha ^ 1 and [SII]AA6716+6731/Ha - 0.9. 
Upon moving away from the nucleus, two clear trends 
in excitation are evident: (1) towards the south- 
east and north-west, i.e., in a direction perpendicu- 
lar to the band, the [OIIIJA5007/H/3 ratio increases 
to around 15 — 20 while [NII]A6583/Ha decreases to 
~ 0.5-0.6, and [SII]AA6716-K6731/Ha also decreases 
to ~ 0.5; (2) along the band, towards the south-west 
and north-east, the behaviour is somewhat inverted, 
with [OIIIJA5007/H/? remaining at a value ^ 10 and 
[NII]A6583/Hq: increasing to values ;> 1.2. The band 
also coincides with the region of high reddening in- 
ferred from the Ha/H/3 ratio, suggesting that this 
region may correspond to a ring or disk of dust-rich 
material. The one-dimensional composite spectra dis- 
cussed in the following sub-section allow a more de- 
tailed investigation of the physical properties of vari- 
ous regions in the EELR. 

3.3.2. Composite Spectra 

Low-dispersion composite spectra are presented in 
Figure ^ for the entire central emission-line region, the 
nucleus (3" x 3"), and the summed emission from the 
three tail regions A, B and C. The wavelength scales 
of the spectra have been corrected to rest wavelength, 
using the mean redshift of the brightest emission lines. 

In Figure ^, we present the summed spectra for the 
Ha object to the north of PKS 2356—61. Its spectrum 
is typical of a star-forming region, with a very young 
stellar population and strong Balmer and [OIIJA3727 
emission, with very little [OIIIJA5007. It is interest- 
ing to note that the two extensions both show some 
emission-line features, although their stellar continua 
are very different. In particular, the southern exten- 
sion has a spectrum very similar to the nucleus of 
the object, while the northern extension has a much 
larger Balmer line equivalent width, displaying very 
little continuum emission. This may indicate that the 
northern extension consists almost exclusively of star- 
forming regions, while the stronger continuum in the 
nucleus and southern extension reveals the additional 



presence of somewhat older stars; a more definitive 
interpretation is precluded by the relatively low S/N. 

Line fluxes were measured using the IRAF routine 
"specfit" ; since these spectra are of low dispersion, 
each profile can be adequately described by a single 
Gaussian component. In clouds B and C, the [SII] 
lines were not fitted due to the noise contribution from 
skylines at 7300A. Table || lists the measured fiuxes 
and l-fj rms errors for all the lines fitted. The errors 
are the formal 1-cr errors obtained from the fits, but 
they do not represent other effects such as sky back- 
ground subtraction uncertainties. Thus their primary 
use is simply to give an indication of the relative ac- 
curacy of the various measured fluxes, particularly 
relevant to the fainter regions. 

In Table ^ we present a number of physically im- 
portant line ratios, corrected for reddening effects us- 
ing several standard assumptions. The intrinsic ratio 
I\-^/Ix2 of the intensities of two lines at wavelengths 
Ai,A2, given an observed ratio I\^ o/I\2 o ^-^d an ex- 
tinction function /(A) (e.g.. Savage & Mathis 1979), 
is calculated following Osterbrock (1989): 

hi_ ^ hi^ gC[/(Ai)-/(A2)] 

where the optical depth at A is ta = C /(A), parame- 
terized by the coefficient of extinction C. Extinction 
due to our galaxy is assumed negligible as a result 
of the relatively high galactic latitude of the source 
(—55° S), thus all reddening is considered to be inter- 
nal to the source. Case B recombination is assumed, 
with the canonical IIa/H/3 ratio of 2.87. Although 
some authors claim evidence of a slightly higher in- 
trinsic ratio of ~ 3.0 — 3.1 when a hard ionizing spec- 
trum is assumed (Halpern & Steiner 1983; Ferland 
& Osterbrock 1985; Tsvetanov & Yancoulova 1989), 
we find that this difference has a negligible effect on 
the present results. Underlying Balmer-line absorp- 
tion is also taken to be negligible, as we discussed 
previously based on the Ha/H/? ratio and the contin- 
uum flux. Finally, it is assumed that reddening for 
all emission-lines can be described by the observed 
Balmer decrement, and that the electron tempera- 
ture Te derived from the [OIII]AA5007/4363 ratio can 
be applied to the Balmer-line emitting region. In 
practice, different ionic species are probably not pre- 
cisely co-incident; however, dust variations within the 
EELR are likely to be less significant than foreground 
reddening, and furthermore the reddening values ob- 
served in PKS 2356—61 are sufficiently small that 



10 



8x10-'= 



6x10- 



ap 4x10-15 

01 

X 
3 

E 2x10-" 



3x10-15 



" 2x10-15 



§ 10-15 




•x; 4x10-1= 

e 

o 

M 

^ 2xl0-i« 

X 
3 



(c) 



-l i I I I l- 

PKS 2356-61: Tail Region 



4000 4500 5000 5500 

Rest-frame "wavelength (A) 



6000 



6500 



Fig. 5. — Low-dispersion spectra (rest wavelength) for: (a) the central region, (&) the nuclear region (3" x 3"), 
and (c) the entire tail region. These spectra were produced by summing all the individual spectra corresponding to 
pixels that have a detectable [OIII] flux in the emission- line image, and thus represent high S/N integrated spectra 
for each region. 
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Fig. 6. — Low-dispersion spectra (rest wavelength) for the faint Ha object to the north of PKS 2356—61. These 
spectra were produced in the same way as those presented in Figure ^, by summing all the individual spectra 
corresponding to pixels with a detectable Ha flux. The top panel (a) shows the total summed spectrum for all 
pixels covering the entire Ha object, including the northern and southern extensions and its central region. The 
second and third panels (b) and (c) display the spectra for the southern and northern extensions, and are shown 
here in order to afford a comparison between the two extensions and the total spectrum, the differences possibly 
being a result of different stellar populations in each region. 
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Table 3 

Line fluxes and errors for the low-dispersion spectra of PKS 2356-61. 



Emission 
line 


Centre 




Tail (Tot.) 


Tail {A) 


Tail (B) 


Tail (C) 


Flux 


1- 


-a 


Flux 


1-0- 


Flux 


1- 


-a 


Flux 


1- 


■a 


Flux 


1-0- 


[OII]A3727 


226.0 


3 


A 


10.9 


0.6 


5.14 


0, 


.87 


3.26 


0, 


.82 


1, 


.03 


0.65 


[NeIII]A3869 


54.3 


2, 


.1 


4.53 


0.51 


2.20 


0, 


.64 














H7 


35.8 


1, 


.7 


2.45 


0.45 


0.47 


0. 


.36 














[OIII]A4363 


20.7 


1, 


.6 


1.88 


0.40 




















HeIIA4686 


19.4 


1, 


.2 


1.84 


0.33 




















H/3 


53.0 


1, 


,2 


2.35 


0.43 


1.33 


0, 


.26 


0.60 


0, 


.29 


0, 


.47 


0.21 


[OIII]A5007 


662.0 


1, 


.8 


38.3 


0.7 


19.4 


0, 


.5 


11.0 


0, 


.4 


4, 


.48 


0.28 


[NI]A5200 


9.95 


1, 


.06 


0.67 


0.41 


0.35 


0, 


.23 














[OI]A6300 


32.7 


1, 


,2 


0.60 


0.35 


0.33 


0, 


.21 














Ha 


229.0 


1, 


.4 


10.1 


0.6 


4.98 


0, 


.26 


2.54 


0, 


.28 


1, 


.68 


0.23 


[NII]A6583 


203.0 


3, 


.7 


7.12 


0.57 


3.02 


0, 


.27 


2.18 


0, 


.30 


1, 


.11 


0.22 


[SII]A6716 


95.0 


2, 


.2 


1.27 


0.63 


0.93 


0, 


.42 














[SII]A6731 


70.3 


1, 


.7 


1.67 


0.58 


0.77 


0. 


.37 















Note. — All values are in units of 1 x 10 ergs ^ cm 



Table 4 

De-reddened line ratios for PKS 2356-61. 



Line Ratio 


Centre 


Tail (Tot.) 


Tail {A) 


Tail (B) 


Tail ((7) 


C 


1.169 


1.189 


0.761 


1.125 


0.647 


HeIIA4686/H/3 


0.39 lil 


0.84 HI 








[OIII]A5007/H/3 


11.9 


16 }I 


14 1^ 
12 


18 ?^ 


9.5 If f 


[OIII]AA5007+4959/4363 


34.5 


22 ?| 








[Olll]A5007/[OI]A6300 


27.4 t! 


87 iJI 


79 146 
<Z 47 






[OIII]A5007/[OII]A3727 


1.90 


2.3 2^2 


2.8 li 


2.2 ?J 


3.4 li 


[OIII]A5007/[NII]A6583 


4.69 ill 


7.9 ?:i 


8.2 n 


7.2 li 


4.9 t:I 


[OIII]A5007/Ha 


A 1 K 4.19 
4.10 4^2 


5.6 li 


4.9 t ? 


6.2 i:i 


q Q 3.6 
3.0 


[OI]A6300/Ha 


1 

U.iO 0^4 


0.06 I j^ 


0.07 III 






[SII]AA6716+6731/Hq! 


0.70 


0.28 I f, 


0.34 lil 






[NII]A6583/Ha 


0.89 li\ 


71 0-79 
^■'^ 0.63 


0.60 HI 


0.86 


0.66 l-ll 


[SII]AA6716/6731 




n 7fi 1-46 
^■'^ 0.36 


J-.^ 0.6 







Note. — The ratio upper and lower limits are derived from the 1-a flux errors. 
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such differences are likely to be undetectable. 

3.4. Physical Conditions in the Gas 

3.4.1. Electron Density and Temperature 

We use the line flux ratios of [SII]AA6716/6731 and 
[OIIIJAA5007/4363 to constrain the electron density 
rie and temperature Tg of the gas. For the central 
region, we find Tf. = 2.2 ± 0.1 x lO^'K and constrain 
He ^ "2,2 — 140 cm"^, with a most probable value of 
He = 80cm~'^. The three tail regions are treated as 
a single region, distinct from the central region, since 
the differences between the tail and central regions 
are generally much greater than within the tail, and 
moreover S/N limitations allow only one value to be 
obtained for the entire tail region. For the tail region, 
we measure Te = 3.2 ±0.6 x lO"* K and constrain Uf, to 
lie in the range ~ 10 — 10^ cm~'^ (with a most prob- 
able value of Uf. ^ 100 cm~'^). The upper and lower 
bounds given for each quantity are derived from the 
errors in the measured ratios in Table ^. We note 
that the values derived from these ratios may not be 
representative of the line-emitting gas as a whole, as 
its inferred ionization structure is model-dependent. 
Therefore, the derived densities/temperatures are in- 
dicative of mean values, possibly integrated over sev- 
eral zones of different ionization and density. 

The [OIII] temperature is comparatively large for 
both the central region and the tail. It is probable 
that the errors on the measurement for the tail region 
are larger than the formal 1-ct errors yielded by the 
flux fitting, although we find that the maximum pos- 
sible uncertainty introduced by the underlying con- 
tinuum results in a lower limit of Tg - 2.1 x lO'^K 
for the tail region. It should be noted that substan- 
tial [OIIIJA4363 emission is observed across the entire 
central region (~ 1 2" x 9" ) ; less than 30% of the flux is 
produced by the inner 3" x 3" centred on the nucleus. 
Thus, the [OIII]AA4363/5007 ratio indicates an elec- 
tron temperature Tg ^ 2 x lO"' K across an extended 
region of > 20kpc; we discuss this interesting result 
in more detail later in this paper. 

3.4.2. Ionization rate 

We derive the ionization rate in the gas, A'ion, using 
the observed Ha luminosity in various regions of the 
EELR, and following standard assumptions such as 
Case B ionization equilibrium conditions (e.g., Oster- 
brock 1989). The various values of Mon are presented 
in Table || (calculated from the Ha fluxes in Table |^). 



The volume of each region is derived from its pro- 
jected spatial extent by assuming an approximately 
ellipsoidal three-dimensional geometry, with a line- 
of-sight depth equal to the mean of its minimum and 
maximum extent on the sky (in fact, regions B and C 
are represented by spherical volumes of diameter d). 
We also tabulate the derived quantity rie/^/^, where 
/ represents the volume filling factor of line-emitting 
gas. It should be noted that this is only an approxi- 
mate indication of the electron density, applicable to 
material in ionization equilibrium; the occurrence of 
shocks can produce post-shock density enhancements 
two to three orders of magnitude higher than the pre- 
cursor gas. 

Later in this paper, we use the values of A^ion to 
investigate various possible ionization mechanisms for 
the tail clouds, as well as the covering fractions of the 
clouds to ionizing radiation, and their geometrical lo- 
cation with respect to the nucleus and the plane of 
the sky. The limits on the electron density from the 
quantity n^f provide valuable additional constraints 
to those derived previously from the [SII] line ratios: 
if shocks are indeed a dominant ionization mechanism 
in the gas, it is likely that the densities suggested by 
the [SII] ratios apply primarily to the post-shock ma- 
terial, while the precursor densities may be closer to 
the lower limits tabulated in Table |[ In the follow- 
ing section we discuss in more detail the question of 
the relative contributions of precursor and post-shock 
emission, through a detailed examination of observed 
line ratios, fluxes and velocity widths, and a compar- 
ison with results from auto-ionizing shock models. 

4. Emission-line Energetics and Excitation 

Here we investigate the applicability of the auto- 
ionizing shock model to the EELR, by examining 
three complementary observable quantities: (a) veloc- 
ity width along the line-of-sight; (b) surface brightness 
for a region of particular velocity width; (c) line ratios 
for a region of particular velocity width. If shocks are 
the principal ionization mechanism in the gas then 
these three observables should be inter-related; this 
allows good constraints to be placed upon the con- 
sistency of the model. The comparison between the 
shock model and the observations is carried out in 
three different ways: 

1. Energetic feasibility. Two aspects are con- 
sidered: (a) The variation of observed luminosity as a 
function of velocity width, for all pixels, as compared 
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Table 5 

Ionization rates in the EELR. 





Extent 


Volume 


A^ion 






(arcsec) 


(kpc^) 




(cm ^) 


Centre 


10 X 8 


2425 


3.14 X 10^^ 


0.176 


Tail A 


6x9 


1365 


6.83 X 10^1 


0.034 


Tail B 




1155 


3.48 X 10""^ 


0.027 


Tail C 




730 


2.30 X lO'^^i 


0.028 



with the relationships expected from the shock model; 
(b) Can the EELR luminosity be accounted for ener- 
getically by shocks with velocities corresponding to 
the observed velocity widths? 

2. Excitation: Dependence on Ad. This com- 
parison is made by examining the excitation state of 
each pixel (quantified using a ratio of two emission 
lines) as a function of its velocity width. Resulting 
trends are compared with those expected on the basis 
of the shock model. This comparison uses the values 
measured for each pixel and is therefore limited to the 
brightest lines, observable across the entire EELR. 

3. Excitation: Diagnostic line ratios. The 

degree of excitation can also be quantified by plot- 
ting pairs of line ratios against each other, specifi- 
cally aimed at separating regions where different ex- 
citation mechanisms operate. In this comparison, in- 
tegrated line fiuxes are obtained by summing spectra 
from different regions of the EELR, in order to obtain 
good signal-to-noise line fluxes for the fainter emission 
lines, some of which are crucial for determination of 
physical properties of the gas. The location of the re- 
sulting points, relative to the loci from shock models, 
yields a powerful means of verifying (independently of 
the previous two comparisons) whether the observed 
excitation state is consistent with the observed veloc- 
ity width if the shock model applies. 

In these comparisons, we assume that the observed 
line-of-sight velocity width Av represents one or more 
shocks with that velocity, rather than a superposition 
of kinematically unrelated systems. Observationally, 
we define Av as the largest line-of-sight velocity dif- 
ference for emission in a given pixel. For strongly 
blended high S/N profiles, we find that the second 
moment of the profile provides a robust means of 



quantifying its velocity width; otherwise if multiple 
peaks are clearly apparent then Av is taken to be the 
largest velocity difference between separate compo- 
nents. The resulting values of Av were "de-blended" 
by the velocity resolution of 70kms^^; however this 
is not a substantial correction. 

4.1. Auto- Ionizing Shock Models 

In order to investigate the viability of auto-ionizing 
shocks as a possible ionization mechanism for the gas, 
Koekemoer (1996) has computed a grid of models ex- 
tending across two of the primary shock parameters: 
the shock velocity and the magnetic field parameter 
B/n^/^, where B is the value of the magnetic field 
strength transverse to the shock, and n is the pre- 
cursor gas density. The gas is assumed to be in the 
low-density limit in all cases. The shock models were 
computed using the non-equilibrium ionization code 
MAPPINGS II, described in Sutherland & Dopita 
(1993), together with the same atomic abundances. 
These abundances are generally the same as those 
in the Solar neighbourhood; the question of possible 
abundance variations is explored further in some of 
the subsequent sections. 

In Figures ^ - |l^ we present results from the shock 
ionization models together with the observational 
determinations of emission-line fluxes and velocity 
widths. The model results plotted are: pure shocks 
(ranging from 150 — 550 km s^^), pure precursor spec- 
tra produced by photoionization from shocks with ve- 
locities in the range 200 — 550 km s~^, and a set of loci 
for the combined spectra of photon-bounded precur- 
sors and shocks in the velocity range 200—550 kms~^. 
Results are also plotted for different values of the 
magnetic field parameter: B/n^^'^ — 0, 1,4/LtGcm^/^. 

An increase in the magnetic parameter reduces 



15 



the compression factor behind the shock, thereby 
increasing the effective ionization parameter in the 
photoionization-recombination region and changing 
the hne ratios accordingly. However, the ionizing 
spectrum of the shock, hence also the precursor emis- 
sion, is insensitive to the value of B/n^/"^. Thus each 
shock velocity produces a single photoionized precur- 
sor spectrum, regardless of different values of the mag- 
netic field parameter. 

4.2. Energetic Feasibility 

We first investigate whether the emission-line re- 
gions may feasibly be powered by shocks with ve- 
locities corresponding to the observed values of At). 
We compare the observations with the modelling re- 
sults by assuming that for a given detector pixel, the 
observed velocity width corresponds to one or more 
shocks along the line-of-sight having only that veloc- 
ity. This is physically justified by the strong depen- 
dence of shock luminosity upon velocity, such that 
the fastest shock in any region can be expected to 
dominate the spectrum. 

The total surface area of shocks along the line-of- 
sight corresponding to a single pixel is parameterized 
by the quantity normalised to the pixel surface 
area. This is necessary because the resolution scale 
is of order 1 kpc, which may be substantially larger 
than the scale of shocks in the gas. Therefore it is 
not possible to directly constrain orientation effects 
for individual shocks, or the possible presence of mul- 
tiple shocks of different velocities along the line-of- 
sight sampled by a single pixel; the inclusion of such 
effects would introduce too many unknown parame- 
ters into the problem. Rather, the approach taken 
here is to examine the degree of consistency of shock 
models with the observations using a minimum of free 
parameters. 

Given the above scenario, the observed flux at each 
pixel is converted to an intensity, in units of ergs~^ 
emitted per cm^ of material along the line-of-sight, 
thus allowing a direct comparison with the intensities 
predicted by shock models. The [OIII]A5007 and H/3 
intensities for all pixels across the EELR are plotted in 
Figure ^ as a function of the observed velocity width 
of each pixel. The models on the plots represent the 
intensity that a pixel of a given velocity width would 
have, given the above assumptions. The models are 
calculated for a precursor density n = 1 cm~^ and 
/a = 1; it should be noted that the predicted intensity 
is linearly proportional to density. 



The models are not intended to be a "fit" to the 
data points. Rather, the primary result from this 
comparison is: 

For any given point (i.e., pixel) in the gas 
with a particular velocity width, shocks of 
a corresponding velocity can account for 
all the flux observed in that pixel for values 
ofnfA^O.1-10. 

There is no reason a priori to expect the shock model 
predictions to be at all comparable to the observed 
intensities. The fact that agreement to better than 
an order of magnitude is obtained for the EELR en- 
ergetics, using physically reasonable approximations, 
provides strong support for the argument that shocks, 
if present, can readily dominate the EELR ionization 
energetics. 

It is clear that the upper envelope of pixel inten- 
sities increases strongly as a function of Ai>, rising 
by an order of magnitude across a factor of two in 
velocity width. In particular, none of the brightest 
pixels display low velocity widths, and this cannot be 
explained through the presence of any kind of observa- 
tional selection effect. Indeed, observational selection 
effects should be more prominent at low intensity lev- 
els: for example, weak lines should be preferentially 
detected at small velocity widths. The fact that all 
the low-intensity pixels display the full range of ve- 
locity widths indicates that these observations are of 
sufficient S/N not to suffer from such selection effects. 
The plots show that the EELR contains pixels with 
(projected) shock velocities continuously ranging be- 
tween ~ 200 and 350 km s~^, while the spread in pixel 
intensities of similar velocity width may indicate dif- 
ferent amounts of line-emitting gas. 

It is also evident that for all values of A-y above 
~ 200kms"\ the observed [OIII]A5007 fiux, relative 
to H/3, remains about a factor of 2 higher than pre- 
dicted by models involving combined shock and pre- 
cursor emission. We show in the following sub-section 
that such a constant offset suggests geometrical ori- 
entation effects, in which an inclination of the shock 
plane to the line-of-sight leads to an observed velocity 
that is lower than the intrinsic shock velocity. 

The energetics alone cannot provide further sig- 
nificant insights into the physical processes along a 
given line-of-sight, since both the shock surface area 
and the density are free parameters. In particular, 
it is difficult to constrain the relative importance of 
emission from the shocked gas versus that of mate- 
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Fig. 7. — Variation of the intensity of (a) H/3 and (b) [OIII]A5007 for each pixel in the EELR of PKS 2356-61, 
as a function of the hne-of-sight velocity width of each pixel. Points from the "Central Region" are indicated by 
filled dots, while those from the Tail regions A, B and Care marked accordingly. Auto- ionizing model curves are 
plotted for "pure" shocks ("Sh"), "pure" precursors, and combined shock and precursor emission ("Sh-|-P"). Note 
that the value of the 5-parameter significantly affects emission from the shocked region, but not from the precursor 
material. 
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rial photoionizcd by the shocks. In the following sub- 
section we investigate line ratio variations, which offer 
a powerful complementary means of constraining the 
relevant physics. 

4.3. Excitation and Velocity Width 

The free parameters of shock surface area and gas 
density can be removed by investigating the varia- 
tion of line ratios as a function of velocity width. In 
particular, a high-excitation line such as [OIIIJA5007 
compared to H/3 provides a good diagnostic of the 
relative contributions of shock and precursor emis- 
sion. In Figure^ we present the ratio [OIIIJA5007/H/3 
as a function of line-of-sight velocity width, for all 
pixels with detected H/? fluxes > 3cr. The values of 
the [OIIIJA5007/H/3 ratio in Figure ^ range predomi- 
nantly between 10 and 20, at all the observed velocity 
widths. However, the points separate clearly into the 
following two regimes, depending upon which part of 
the EELR they represent. 

Firstly, the central region and the tail region B 
both display velocities up to ~ 350kms~^ (typically 
in the range 250—300 km s~^); their observed values of 
[OIII]A5007/H/3 correspond to a combination of emis- 
sion from photoionizcd precursor material (kinemati- 
cally associated with shocks) and cooling material be- 
hind shocks of velocities ~ 400 — 550 km s^^. These 
shock velocities become consistent with the observed 
velocity widths upon invoking the standard projec- 
tion factor of 1/ ^3 for an ensemble of line-of-sight 
shock inclination angles (thereby producing observed 
velocities only ~ 50 — 60% of the intrinsic shock ve- 
locities). 

Secondly, the tail regions A and C display ve- 
locity widths almost exclusively <^ 200kms~^ and 
predominantly in the range 100 — 150kms~^, much 
too small to be kinematically associated with shocks 
of the velocities required to produce the observed 
[OIII]A5007/II/3 ratios (unless highly improbable geo- 
metric situations are invoked). Thus, in the context of 
auto-ionizing shocks these regions correspond to gas 
photoionizcd by shocks, but not kinematically associ- 
ated with such shocks. 

If the line-emitting gas in the central region is re- 
lated to infalling material, for example gas captured 
tidally through an interaction with a companion ga- 
laxy, then such gas clouds are likely to virialize on 
significantly shorter timescales than the overall life- 
time of the source. If shocks in this region are the re- 



sult of collisions between individual clouds, then the 
virialized velocity distribution must be deprojected 
to correspond to the observed velocities. We use the 
standard deprojection factor of 1/^/3 for an ensemble 
of clouds that possess random velocity vectors. The 
actual situation is likely to be more complex: all the 
clouds would not be moving at identical speeds; shock 
velocities would depend upon the relative trajectories 
of the clouds involved; and the cloud velocity vectors 
are not likely to be completely randomized. However, 
it is clear that agreement with the observations is 
substantially improved if the observed velocities cor- 
respond to ~ 50 — 60% of their intrinsic values, not 
unreasonable for an ensemble of clouds in the process 
of virialization. 

4.4. Diagnostic Line Ratios 

The final comparison between observed emission- 
line properties and those predicted from shock mod- 
els, involves the use of "diagnostic" line-ratio plots 
(Baldwin et al. 1981; Veilleux & Osterbrock 1987; 
Dopita & Sutherland 1996). Such plots are designed 
to allow the separation of different excitation mech- 
anisms, thereby providing useful insights into the 
physics of a particular object. In the present pa- 
per they are used to separate pure shocks from pure 
precursors and shocks with "full" (photon-bounded) 
precursors, as well as discriminating between results 
expected for shocks of different velocities and differ- 
ent values of B/n^^^. Hence, these diagrams pro- 
vide an independent means of investigating the con- 
sistency of shocks with the observations. The ratios 
are taken from Table ^ and represent the following 
distinct regions: (1) the Central elongated emission- 
line region; (2) the total Tail region, representing the 
tidal "tail" structure extending away from the central 
region; (3) the three individual "tail" sections A, B 
and C, for which fluxes of some of the brighter lines 
were obtained. 

The relevant ratios and 1-a error limits for these 
regions are presented on the diagnostic line ratio plots 
in Figures ^ - |ll|, together with the loci expected from 
shock models. The 1-a limits are purely indicative of 
the statistical errors of the measured fluxes; they do 
not include effects such as line ratio variations across 
individual regions. Note that fluxes for some of the 
fainter lines could not be obtained in all three tail 
regions A, B and C, thus they are omitted from some 
of the plots. A comparison of the observed data with 
the models allows a number of physical properties to 
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Fig. 8. — Variation of [OIIIJA5007/H/3 as a function of observed velocity width, for all pixels of sufficient S/N 
in the EELR of PKS 2356—61. The data points and theoretical shock loci are represented the same way as in 
Figure ^ Two families of shock loci are plotted (both extending up to intrinsic shock velocities of 550kms~^): one 
for a line-of-sight inclination of 0° (i.e., the shock velocity is identical to the observed line-of-sight velocity width, 
extending up to 550kms~^), while the other corresponds to the velocity widths expected for a random distribution 
of line-of-sight inclinations; (thus the observed ensemble of shock velocities are reduced by the "mean" projection 
factor of l/\/3, and only extend up to 300kms~^). 
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be inferred about the gas and allows a good test of 
the self-consistency of the the shock model, which we 
discuss here for various line ratios. 

[OIII]A5007/H/3 vs. [OII]AA3727/[OIII]A5007 (Fig. 
^). — The relatively high values of the ratio of 
[OIIIJA5007/H/? imply the presence of photoionized 
precursors, since the shocked material alone produces 
comparatively weak [OIIIJA5007 emission. The ob- 
served data points are consistent with combined shock 
and precursor emission, with shock velocities in the 
range ~ 400 — 550kms~^ and a magnetic parame- 
ter value B/v}/"^ ~ 1/iGcm^/^. However, the points 
for A and C lie somewhat closer to "pure" precursor 
emission than the region B and the Central region; 
it should also be noted that uncertainties in the red- 
dening correction could either increase or decrease the 
[OII]A3727/[OIII]A5007 ratio by - 30 - 50%. 

[OIII]A5007/H/3 vs. [NII]A6583/Ha (Fig. |^). — 
The positions of the [NII]A6583/Ha points with re- 
spect to the model curves are remarkably consistent 
with those in Figure even though the ratio of 
[NII]A6583/Hck is observationally independent of the 
[OII]A3727/[OIII]A5007 ratio. This agreement offers 
strong support for the feasibility of the shock model. 
In particular, Figure 1^6 also suggests that the emis- 
sion from the gas is due to shocks combined with 
photon-bounded precursors, and that the magnetic 
field parameter is likely to be B/v}/^ ~ 1/xGcm'^/^ 
rather than zero. This is a typical value of the 
magnetic parameter expected for ISM gas (Dopita & 
Sutherland 1995). It is also interesting to note that 
the agreement with Figure implies that nitrogen 
abundances are not likely to be substantially differ- 
ent from Solar values (i.e., within a factor 2 — 3). 

[OIII]A5007/H/3 vs. [SII]AA6717-h6731/Ha (Fig. 
^0|a). — The tail regions B and C are absent from 
this figure, since the [SII] emission from these re- 
gions is of too low S/N to be useful. Thus the only 
fluxes plotted here are from the Central region, from 
the tail region j4, and from the summed spectrum of 
all three tail regions. The results apparent on this 
diagram are very similar to those of Figure ^ i.e., 
the ratios are dominated by emission from radiation- 
bounded precursors photoionized by shocks of veloc- 
ities in the range 450 — 550kms~^, with a magnetic 
parameter B/ri^/'^ ^ l/iGcm"^/^. In particular, the 
emission from the tail region A (and the tail as a 
whole) appears dominated by emission from precur- 
sors rather than post-shock material, while the cen- 
tral region is, as before, consistent with post-shock 



emission together with kinematically associated pho- 
toionized precursors. We also note that the relative 
consistency of these results with those of Figure || 
implies that additional effects (such as collisional de- 
excitation) are unimportant, thus the gas is predom- 
inantly in the low-density limit (rig <^ 10'^ cm""^). 

[OIII]A5007/H/3 vs. [OI]A6300/Ha (Fig. |l^6). — 
These ratios are also consistent with ionization in- 
volving shocks of velocities ^ 400kms~^. As in Fig- 
ure the [OIJA6300 emission in the tail is rel- 
atively weak and is primarily from region A, the 
location of which suggests emission dominated by 
"pure" precursor material. The Central region lies 
substantially closer to the loci comprising combined 
emission from precursors and shocks with velocities 
~ 400—450 km s^^, and a weak value for the magnetic 
parameter. It should be noted that despite the im- 
portance of [OIJA6300 as a diagnostic tool for discrim- 
inating between optically thick and thin photoioniza- 
tion conditions, its modelling entails substantial un- 
certainties. For example, the formation of molecules 
in the partially ionized region can lower the electron 
temperature and significantly reduce the [01] emis- 
sion (Dopita & Sutherland 1996), and this is not cur- 
rently accounted for in the MAPPINGS II code. Pre- 
dicted model values of the [OIJA6300 line emission 
tend to vary by a factor ~ 2 — 3 between different 
ionization codes (Ferland et al. 1995), a substantially 
larger variation than for the higher ionization states of 
oxygen. The discrepancy between the MAPPINGS II 
models and the observations is also of this order, thus 
it is possible that incorporating this effect into the 
models would make Figure [lO| b more consistent with 
Figure |lO|a, although such modifications are well be- 
yond the scope of the present work. Despite this un- 
certainty, however, the results from the [OI]A6300/IIq; 
ratio are still remarkably consistent with those ob- 
tained from the other ratios. 



[OIII]A5007/H/3 vs. HeIIA4686/H/3 (Fig. lli) 



The HeIIA4686 fiux for the central region agrees well 
with the scenario inferred from the other ratios, i.e., 
combined emission from shocks of velocity ~ 400 — 
500kms~^ with radiation-bounded precursors, and 
magnetic field parameters of order ^ 1 /iC cw?/"^ . The 
value of IIeIIA4686/H/3 for the tail region is approx- 
imately consistent with pure precursors photoionized 
by shocks of velocity ^ 400 — 450kms~^; the value 
of the ratio is somewhat higher (about a factor of 2) 
than expected from the models, but it should also 
be noted that the observational uncertainty for this 
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Fig. 9. Diagnostic ratios of: (a) [OIII]A5007/H/3 vs. [OII]AA3727+3729/[OIII]A5007; (6) [OIII]A5007/H/3 vs. 
[NII]A6583/Ha. All the model curves extend to 550 kms~^, with velocities plotted at increments of 50 kms~^. Loci 
for pure precursor emission commence at 200kms~^, while those involving shocks commence at 150kms~^. Note 
that the value of the B-parameter significantly affects emission from the shocked region, but not from the precursor 
material. 
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ratio in the tail region is comparatively large. Thus, 
these results are overall approximately consistent with 
those obtained from Figures ^ and 

[OIII]A5007/H/3 vs. [OIII]A4363/[OIII]A5007 (Fig. 
^b). — Upon comparison with the data presented in 
Dopita & Sutherland (1995) it is clear that the elec- 
tron temperature for the gas in this galaxy is among 
the highest seen in radio galaxy emission line regions. 
It is likely that uncertainties in the underlying con- 
tinuum features have led to an overestimation of the 
[OIIIJA4363 flux; however, the largest possible error 
introduced in this way corresponds to an overestima- 
tion of the [OIII]A4363 flux by a factor of - 2-3. Re- 
ducing the observed [OIIIJA4363 flux by this amount 
would still provide good agreement, not only with the 
fluxes predicted by the shock model, but also with the 
particular scenario inferred from the other diagrams 
— combined precursor and post-shock emission from 
shocks of velocities 400 — 550 km s^^. By compar- 
ison, nuclear photoionization models typically under- 
predict this ratio by factors > 50—100 (e.g., Stasihska 
1984; Tadhunter, Robinson & Morganti 1989). Such 
discrepancies can be reduced to some extent through 
the use of lower than solar abundances (Ferland & 
Mushotzky 1982; Ferland & Netzer 1983) or the in- 
troduction of matter-bounded clouds having various 
distributions of thickness (Morganti et al. 1991; Vie- 
gas & Prieto 1992; Binette et al. 1996). However, in 
order to compare those results with the models pre- 
sented here, we would need to introduce similar free 
parameters, such as different metallicities, and vary- 
ing distributions of precursor optical depths or shock 
velocities along the line-of-sight, and such modifica- 
tions are outside the scope of the present work. 

In summary, the following scenario can be inferred 
from the comparison of three observationally inde- 
pendent fine ratios, [OIII]A5007/H/3, [NII]A6583/Ha 
and [OII]A3727/[OIII]A5007: the ratios in the cen- 
tral region and the tail region B are consistent with 
combined emission from shocks with velocities in the 
range 400 — 550 km s^^ and magnetic field parame- 
ters of order B/n^^^ — 1/iGcm'^/^, together with 
radiation-bounded precursors photoionized by these 
shocks. The ratios in the tail clouds A and C tend to 
be more indicative of emission dominated by precur- 
sors, photoionized by shocks with velocities of order 
400 — 550kms~^. The fact that consistency between 
the observations and the models can be achieved 
in these ratios, without invoking additional free pa- 
rameters such as "partial" precursors (i.e., matter- 



bounded), density differences, or abundance varia- 
tions, is a very strong point in support of the viability 
of auto-ionizing shocks as an ionization mechanism for 
the EELR. 

The examination of two additional line ratios, 
[OI]A6300/Ha and [SII]AA6717+6731/Ha, suggests 
that the temperature in the extended weakly ionized 
region may be somewhat lower than expected, possi- 
bly as a result of the formation of molecular species. 
This is the only modification required that is not pos- 
sible to incorporate into the current model. It is also 
apparent that the gas is predominantly in the low- 
density limit. Finally, the shock model accounts for 
the ratio [OIII]A4363/[OIII]A5007 to within a factor 
of 2 — 3 (or closer, if uncertainties due to the underly- 
ing continuum are taken into consideration) , and this 
is more than an order of magnitude superior to the 
estimates generally obtained from nuclear photoion- 
ization models. 

5. Discussion 

We now consider the feasibility of different possible 
physical scenarios inferred from the observations, in 
terms of the spatial geometry, dynamics and ionizing 
mechanisms operating in the EELR. For simplicity, 
we discuss separately the two most plausible sources 
of ionizing photons: the AGN or shocks. While it is 
possible that both of these mechanisms are operat- 
ing to some extent, perhaps even to different degrees 
in different parts of the EELR, in the present work 
our aim is rather to investigate the dominant mecha- 
nism responsible for the bulk of the EELR energetics 
across the largest spatial scales. This greatly reduces 
the number of possible combinations of free param- 
eters, thereby allowing more useful constraints to be 
obtained on the physics of the line-emitting gas. 

5.1. Structure and Geometry of the Tail EELR 

Firstly, the structure and physical properties of the 
three clouds in the tail region are investigated using 
the ionization rates presented in Table ^, and the ob- 
served projected spatial distribution of the clouds on 
the sky. This enables constraints to be placed on their 
intrinsic distances from the nucleus, together with re- 
lated physical properties such as the covering fraction 
and volume filling factor of line-emitting gas in these 
clouds. 

We consider a scenario in which the three tail 
clouds are primarily photoionized by radiation escap- 
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ing from the central region. This is motivated by 
the relatively low linewidths in the tail clouds: for 
example, the velocity width across all of cloud C is 
< 150 km s^^, much too small to be directly associ- 
ated with the shock velocities required to ionize the 
gas (unless all the shocks are highly inclined to the 
line of sight, which is very improbable). We later dis- 
cuss the possible effect of supplementary ionization, if 
shocks with velocities ~ 200 — 300 km s~^ are present 
within the tail clouds A or B. 

If the primary source of ionizing photons is con- 
sidered to lie towards the centre of the galaxy, this 
has the additional advantage that the derived tail 
cloud properties are relatively insensitive to whether 
the radiation is produced by the AGN or a collection 
of shocks throughout the central region. We note, 
however, that the inferred properties of the central 
EELR depend strongly upon the scenario invoked for 
production of the ionizing radiation, and in the next 
sub-section this is used to investigate the relative vi- 
ability of shocks or the AGN as the primary source of 
ionizing photons. 

We denote by A^tot the total rate of ionizing pho- 
tons produced by the radiation source; a fraction Tccn 
of these photons are absorbed by the entire central 
EELR, while the remaining flux of escaping ionizing 
photons is Nf.sc = (1 ~ -^^con) A^tot- For each tail cloud 
we use the values of presented in Table ^ to con- 
strain A^csc- The covering fraction of material in each 
tail cloud is represented by T , where T — \ indicates 
total absorption of all incident ionizing radiation by 
the cloud. The observed projected separation A be- 
tween each tail cloud and the nucleus is related to 
their intrinsic separation R by the angle i, defined by 
i? cx A/ cosi (thus ? = 0° corresponds to a geometry 
where the line between the nucleus and the tail cloud 
lies in the plane of the sky). 

In Table || we present N'^^^c = A^esc-^cos^ i, which 
for each tail cloud corresponds to the inferred total 
rate of ionizing photons escaping from the ionizing 
source, required to produce the observed values of 
A^ion (taking into account the parameters !F and i for 
each cloud). The projected distance between each 
cloud and the nucleus, together with the estimated 
spatial extent of each cloud, are used to determine 
the solid angle subtended by each cloud at the nu- 
cleus. We also define a quantity R' , representing the 
distance between the nucleus and tail clouds A and B 
that would be required to produce the same value of 



N'^sc calculated for tail cloud C: 

Thus, R' parameterizes the intrinsic distance R be- 
tween a tail cloud and the nucleus, together with the 
covering fraction !F of the cloud, relative to Tc (the 
covering fraction of cloud C) and ic, the angle be- 
tween the vector from the nucleus to cloud C and the 
plane of the sky. Cloud C is chosen as a reference 
because it is the furthest from the nucleus, thereby 
minimizing the effects of uncertainties in estimates of 
its spatial extent and its separation from the nucleus. 
Although some uncertainties are involved in these 
derivations, for example approximations concerning 
the precise three-dimensional structure of each tail 
cloud, it is still possible to investigate the contribu- 
tion of effects such as the inclination and covering 
fraction on our estimates of Aosc, under the following 
three physically different scenarios. 

Firstly, if the covering fractions J- are approxi- 
mately the same for each tail cloud, then the values 
of N'f^sc and R' suggest a geometry where tail cloud C 
is closest to the plane of the sky (with respect to the 
nucleus) and tail cloud A is furthest out of the plane 
of the sky. More specifically, if the vector between 
the nucleus and tail cloud C lies approximately in the 
plane of the sky (i <^ 20°), then i ^ 50° for the vector 
between the nucleus and tail cloud B, and i ~ 65° for 
the vector to cloud A. This would imply that all three 
tail clouds are at almost the same intrinsic distance 
from the nucleus {R ~ 25kpc), moving on an orbit 
with low ellipticity. Interestingly, this could be con- 
sistent with the observed kinematics: cloud C shows 
the greatest velocity offset with respect to the galaxy 
nucleus, cloud A displays a velocity almost identical 
to systemic (as would be expected if its tangential or- 
bital velocity is strongly projected), and cloud B dis- 
plays an intermediate velocity. One weakness of this 
scenario is that it does not readily account for the ob- 
served velocity split in tail cloud B\ we later discuss 
the kinematics in more detail and consider possible 
extensions to this scheme. 

A second possibility is that the covering fraction of 
the tail clouds increases as a function of radius. If all 
three tail clouds have the same value of i, then the 
covering fractions for clouds B and C relative to cloud 
A are Tb I.STa and Tc S.SJT^ respectively. 
This might possibly be accounted for in terms of ra- 
dial differences in the external medium around each 
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Table 6 

Flux of Escaping Ionizing Radiation. 





Distance A 




R' 




(arcsec) 




(kpc) 


Tail A 


6" 


5.8 X 10''^2 


22.2 


Tail B 


10" 


1.1 X 10''^3 


26.5 


Tail C 


15" 


2.3 X 10^3 


27.9 



cloud (as discussed in Binette et al. 1996, or Prieto 
et al. 1993). For example, if the covering fraction is 
directly proportional to the volume filling factor, and 
if this depends in turn upon the external pressure, 
then a radial decrease of pressure as ~ r~^'^ may pro- 
duce the observed change. However, this is inconsis- 
tent with the radial profiles of or steeper generally 
obtained for ellipticals and similar systems (Kochanek 
1994b; de Zeeuw & Franx 1991 and other references 
therein). Furthermore, there exist no detailed phys- 
ical models describing such effects produced by the 
surrounding medium upon the volume filling factor 
of material within a line-emitting region, hence this 
scenario is difficult to quantify further. 

Finally, in the context of photoionization by the 
AGN, the differences in A^osc might be explained 
by different amounts of radiation reaching the three 
clouds, as a result of anisotropic emission. For exam- 
ple, cloud A may lie further toward the edge of an 
ionizing "cone" , while clouds B and C may be closer 
to its axis, thus more of the ionizing radiation from 
the AGN is obscured towards cloud A compared with 
B or C. Unfortunately, in this scenario the possible 
combinations of cloud positions and cone opening an- 
gles and relative degrees of obscuration are almost 
entirely arbitrary — the observed structure may be 
matched with several completely different cloud posi- 
tioning configurations and cone angles, which do not 
readily allow independent verification through other 
observational properties. Perhaps more importantly, 
a highly anisotropic ionization cone is not easily rec- 
onciled with the uniform radial flux distribution of 
the central EELR. 

We adopt here the simplest parameterization for 
the observed differences in A^esc, assuming that they 
are primarily due to changes in a single variable, ei- 
ther i or J^. If we consider i to be the free vari- 



able, setting to be the same for all three clouds, 
then the tabulated values of R' (or A^^c) imply that 
the clouds are all at approximately the same distance 
from the nucleus {R'/cosic ~ 25 kpc). Hence, the 
clouds would all be in approximately the same envi- 
ronment, thereby providing an independent physical 
justification as to why should be similar for all the 
clouds. The alternative possibility, considering to 
be the dominant free parameter and fixing i to be the 
same for each cloud, is inherently more complex since 
it demands additional assumptions concerning differ- 
ences in the environment around each cloud, as well as 
a detailed description of how changes in the environ- 
ment translate into differences in covering fraction. 

ThTis we consider ilic different values in A'^gf, to 
arise primarily from different values of i, implying 
that the tail clouds are at similar distances from the 
nucleus (small variations may be accounted for by dif- 
ferences in !F, related to the detailed shapes of the 
clouds). This scenario is appealing for several reasons: 
(1) the number of free parameters is minimal; (2) the 
inferred values of i are comparatively insensitive to 
the spatial structure of the source of ionizing photons 
(i.e., AGN or shocks), provided that they originate 
closer to the nucleus than the tail cloud itself; (3) the 
derived tail cloud locations yield a dynamical model 
that can be rigorously tested through a comparison 
with the kinematic structure. 

5.2. Energetics of the Ionizing Source 

The tail cloud fluxes provide constraints upon 
global properties of the central EELR, such as the 
fraction of ionizing photons that are absorbed, given a 
speciflc scenario for the ionization mechanism. In this 
section we examine the energetics required to account 
for the emission from the central EELR and the tail 
clouds, in the context of either shocks or AGN pho- 



26 



toionization. We also extrapolate the ionizing spec- 
tra into the optical regime and compare the expected 
fluxes with the observed low-dispersion spectra. 

Since the tail cloud C is considered to be entirely 
photoionized by external radiation, its value of N'^^^ 
can be used to infer directly the flux of ionizing pho- 
tons, Ncsc{v > 1^0)1 escaping from the ionizing source: 
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2.3 X 10 

J-c cos^ ic 



Here — 13.6 eV is the threshold ionization fre- 
quency for Hydrogen, which we use to obtain the 
total ionizing flux from ionization models (thus al- 
lowing optical fluxes to be computed and compared 
with the observed spectra). We consider the general 
geometry of ionizing photons escaping isotropically 
from the central region, before discussing other pos- 
sibilities such as selective absorption towards certain 
directions (or possible contributions from local shocks 
in the tail clouds). 

5.2.1. Photoionization by the AGN 

We first discuss the EELR energetics in the con- 
text of "unified schemes", in which all the ionizing 
photons are considered to originate from the AGN. 
The total rate of ionizing photons absorbed by the 
central EELR is given by the corresponding value of 
Mon in Table |; an upper limit on the total rate of 
escaping ionizing photons, implied by the tail cloud 
fluxes, is given by the highest value of TV^c (that of 
tail cloud C) in Table Hence we obtain directly an 
upper limit on ^cen, the only unknown parameters 
being the values of T and i for tail cloud C: 
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0.733 



Tc cos^ ic 



If the AGN produces a power-law spectrum of the 
form cx z^" (e.g., Robinson et al. 1987) then the 
value of A^ion for the central EELR, combined with 
the upper limit on Teem yield direct lower limits on 
the expected flux from the AGN in the optical regime. 
For example, setting a = —1.5 implies a flux density 
Fa ^ 5.3 X lO^i^ergs^^cm-^ A-i at A - 4000 A, 
which is at least an order of magnitude higher than 
the observed upper limit on the possible featureless 
component in the nuclear continuum (the absorption 
feature at A = 3935 A in the spectrum of the cen- 
tral 3x3", presented in Figure displays a mini- 
mum fiux level of F\ — 5.1 x lO^^^ergs^"'^ cm^^A^^). 



The most likely interpretation of this deficit involves 
a high degree of obscuration toward the AGN, which 
would be in agreement with the lack of observed nu- 
clear broad emission lines in this object and also gen- 
erally consistent with the "unified" scenario of ob- 
scuration along our line-of-sight towards the nuclei of 
narrow-line radio galaxies. 

The difficulty with accounting for the EELR en- 
ergetics predominantly by radiation from the AGN 
is that the central EELR displays no morphological 
signs of biconical photoionization, but instead has a 
very uniform radial distribution, and is somewhat ex- 
tended along a position angle apparently unassociated 
with the radio axis (offset by ^ 40°). This morphol- 
ogy is also evident in the narrow-band image previ- 
ously published by Danziger & Focardi (1988). Fur- 
thermore, the radial density distribution required to 
account for the tail fluxes is too flat compared with 
physically reasonable models of elliptical galaxies, as 
already discussed in sub-section 5T. Thus, we are 
motivated to investigate another possible mechanism 
for producing the bulk of the ionizing radiation, by 
means of shocks within the EELR. 

5.2.2. Ionization by Shocks 

In this sub-section we discuss the feasibility of ac- 
counting for the EELR energy budget primarily in 
terms of shocks, rather than radiation from the AGN. 
Once again, tail cloud C is considered to be photoion- 
ized, thus we assume that its recombination line lu- 
minosity follows the same scaling relationships as de- 
rived for precursor emission by Dopita & Sutherland 
(1996). 

We consider flrst a geometry where all the ionizing 
radiation is produced by shocks occurring in the cen- 
tral EELR (neglecting for the moment the possibil- 
ity of second-order contributions from lower-velocity 
shocks within the tail clouds A or B ). The physical 
significance of the quantity Tccn in this scenario is 
that it represents the fraction of ionizing photons ab- 
sorbed by precursor material (of density nccn) directly 
associated with the shocks in the central EELR. For 
simplicity, the shocks are all assumed to correspond 
to a single velocity Vg; thus we may directly apply the 
luminosity scaling relationships for shocks and precur- 
sors given in Equations (3.4) and (4.4) respectively of 
Dopita & Sutherland (1996). (together with ratios 
presented in their Tables 8 and 10 to scale the lumi- 
nosities from H/3 to Ha, in order to afford a compar- 
ison with the higher S/N observed Ha fluxes). This 
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yields the expected total observable Ha flux produced 
by the shocks and precursors in the central EELR, 
necessary to produce the required number of escap- 
ing ionizing photons to photoionize the tail cloud C 
(having a density nc)- 

1.68 X 10-1'* 



TcCOS^ ic{l - Tccn) 
0.13 

0.787 



nc 



Vs 



100 kms- 



ergs ^cm ^ 



Although the scenario considered here is compara- 
tively simple, involving shocks of a single velocity, 
and assuming precursors of uniform density, it ac- 
counts remarkably well for the total observed EELR 
luminosities. For example, if Tc cos^ ~ 1 (thus 
-^cen ^ 0-6), then the expected value of Fna.cen, as 
required to photoionize cloud C through photons es- 
caping from shocks, agrees well with the observed 
central EELR Ha flux of 2.29 x lO"*'' ergs"! cm'^ 
(see Table ^) . The agreement can be improved if the 
average density of precursor material in the central 
region is a factor of ^ 2 — 3 times lower than that 
of cloud C, which would physically correspond to the 
scenario of material in the central region being more 
disrupted and diffuse than the cooler gas in the tail 
region. Furthermore, if J-c cos^ ic «C 1 then larger 
density differences can also be accommodated. 

The expected continuum flux density in the opti- 
cal regime from the shocks and precursors in the cen- 
tral region can be obtained directly from the MAP- 
PINGS n model output flies, and compared with the 
observed limits on non-stellar emission. We flnd that 
at A ~ 4000 A, the expected total underlying contin- 
uum emission from shocks and precursors in the cen- 
tral EELR (taking Vs = 500 kms'^ and J^cen ~ 0.6) 
is Fx ^ 1.5 X 10"^'' erg s"^ cm-^A"*; see also Figure 5 
in Sutherland et al. (1993). This is well within the ob- 
served upper limit of i^A = 2 x lO"'^^ erg s""'^ cm-^A"-*^ 
in the absorption minimum at A = 3935 A in the to- 
tal spectrum of the central EELR (Fig. ^) , thus the 
required shock luminosity does not contradict the ob- 
servational limits upon the non-stellar continuum. 

Some shocks may also possibly be present within 
the tail region itself; in particular, the double velocity 
system in tail cloud B is suggestive of two kinemati- 
cally and spatially related emission-line systems. The 
geometry of tail cloud C with respect to B is such that 
it should intercept ~ 15% of any ionizing radiation 
escaping from cloud B; a similar argument to the one 



above reveals that ionizing photons escaping from any 
shocks in cloud B can account only for ^ 10 % of the 
observed emission from cloud C. Thus, the ionizing 
photons from any shocks present within the tail re- 
gion are more likely to be absorbed by their associated 
precursors instead of escaping to ionize other parts of 
the EELR. This would suggest that the volume fllling 
factor of material in the tail clouds is closer to unity 
than in the central EELR, again consistent with a 
scenario in which the central material is strongly dis- 
rupted, while the tail is still relatively quiescent. 

5.3. Shocks: Physical Characteristics 

We have so far demonstrated that the EELR mor- 
phology, excitation, kinematics and luminosities can 
be successfully accounted for if the energy budget is 
dominated by autoionizing shocks, predominantly oc- 
curring within the central EELR, and that this sce- 
nario requires less unknown parameters than the al- 
ternative scheme of photoionization by radiation from 
the AGN. Here we discuss the viability of the shock 
hypothesis in terms of the overall physical properties 
required of the shocks and their precursors. 

The spatial structure of the gas is clearly an im- 
portant consideration in discussing the occurrence of 
shocks. We point out that the volume filling factor / 
of the line-emitting gas, as commonly inferred from 
the density-sensitive [SII] doublet, is only a reliable 
density indicator if the emission is dominated by ma- 
terial with approximately uniform density. Thus, the 
values of rie derived from the [SII] ratio, combined 
with the constraints on rief^^^ derived from the ion- 
ization rates (Table ||) would suggest / ~ 5x 10"^ for 
the emission-line gas in the central region. While this 
is generally consistent with values inferred by other 
authors (e.g., van Breugel et al. 1985, 1986; Heck- 
man et al. 1989), this interpretation is valid only if 
the emission is entirely from material photoionizcd by 
an AGN. 

On the other hand, as we have already alluded 



to in Sub-section 3.4, if the emission is primarily re- 
lated to shocks then the [SII] ratios would be domi- 
nated by emission from the post-shock cooling region, 
where the typical densities are two to three orders of 
magnitude higher than in the precursor zone. For 
example, shocks with Vs ~ 450 — 500 kms^^ and 
B/n^/^ ^ 1 — 2/iGcm^/^ have post-shock regions 
200 — 400 times denser than their precursors, im- 
plying precursor densities riccn <^ 0.1 — 1 cm^^ (if the 
density value for the central region derived in Sub- 
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section 3.4 applies primarily to post-shock material). 
This would suggest values of / within an order of mag- 
nitude from unity for the unshocked precursor mate- 
rial in the central EELR, thereby profoundly altering 
the physical description of the gas. 

The radiative cooling timescale for material in the 
post-shock region is typically ^ 10'* yr, too short for 
a single shock event to account for the ionization of 
the gas during the lifetime of the radio source (gener- 
ally inferred to be ~ 10® — 10^ yr, for powerful ra- 
dio galaxies, e.g., Alexander & Leahy 1987). One 
possibility for continuous deposition of energy into 
shocks is a direct association between the radio jet 
and the line-emitting plasma, as postulated for Cen- 
taurus A (Sutherland et al. 1993), and also more re- 
cently for PKS 2250-41 (Clark et al. 1997). However, 
PKS 2356—61 displays very little association between 
the EELR and the radio axis. In fact, the EELR mor- 
phology and kinematic structure are more consistent 
with accretion onto the host galaxy; thus we propose 
for this object a scenario involving shock formation 
that is driven by collisions between infalling gas, and 
material that is already settling in the host galaxy. 

Using the combined expressions for precursor and 
shock emission as a function of shock velocity, to- 
gether with the observed Ha flux from the central 
EELR, we find that the total surface area of shocks 
in the central EELR required to produce the observed 
emission is: 



A. 



1.44 X 10^ 



0.787 



^cm 



Vs 



Vs 



100 km s- 



100 kms- 



-2.28 



kpc^ 



(where J^cen represents, as before, the fraction of ion- 
izing radiation produced by shocks that is absorbed 
by precursors within the central EELR). Given a pre- 
cursor density neon ~ 0.1 cm^"^, with J-ccn ~ 0.6, and 
shock velocities Vs — 500kms~^, the total mass flux 
through the shocks is M 1.26 Mq kpc~^ yr~^. We 
note the possibility of "reprocessing" of post-shock 
material by further shocks during the lifetime of the 
source, in which case the mass flux through the shocks 
need not necessarily be a direct indicator of the total 
mass of gas in the central EELR. 

5.4. Dynamics and Origin of the Gas 

The strong kinematic and spatial association be- 
tween the clouds in the "tail" EELR suggests that 



they have a common dynamical history; here we dis- 
cuss their origin in the context of gas accretion by the 
host galaxy, as a result of an interaction with a gas- 
rich companion. This scenario is favored on the basis 
of several pieces of observational evidence: (1) the 
presence of the Ha object towards the north of the 
source; (2) the faint bridge of emission-line gas con- 
necting the Ha object with the host galaxy; (3) the 
fact that this faint bridge apparently lies on the same 
trajectory as the "tail" EELR to the south of the host 
galaxy. Furthermore, we rule out the following alter- 
native explanations for the gas origin: (1) ejection of 
gas along the radio axis is inconsistent with the lack 
of any observed association between the EELR and 
the radio axis; (2) condensation from a cooling flow 
in the ICM is not favored because of the compara- 
tively high observed velocities, and also the fact that 
the large-scale EELR is concentrated in a "tail" , not 
a diffuse halo around the host galaxy as expected for 
cooling flows (Heckman et al. 1989; Baum et al. 1992; 
Fabian 1994). 

Numerical studies of gas dynamics during galaxy 
encounters generally show a tendency for the gas dis- 
tribution to bifurcate into two components (Combes 
et al. 1990; Hernquist & Mihos 1995; Barnes & Hern- 
quist 1991, 1996): a compact, centrally condensed re- 
gion, and extended "tidal tails". While the detailed 
morphologies and dynamics clearly depend upon the 
interaction geometry, this general bifurcation is phys- 
ically understood as a direct consequence of the de- 
creased dynamical timescale towards the centre of the 
accreting galaxy, the gas having been driven into the 
central region through gravitational torquing in re- 
sponse to the tidal forces (Binney & Tremaine 1987; 
Barnes & Hernquist 1992): material within the in- 
ner few kpc experiences several orbital crossings and 
begins to virialize while material extended on scales 
^ 10 kpc is still undergoing infall or large-scale orbital 
motions. 

The "settling" of gas towards the centres of el- 
liptical galaxies has also been alluded to in several 
observational studies (e.g., Bertola & Bettoni 1988; 
MoUenhoff & Bender 1990), together with more de- 
tailed simulations of the kinematic evolution of the 
gas (Habe & Ikeuchi 1985, 1988; Varnas et al. 1987). 
A typical feature of such simulations is the presence 
of large velocity shears in the gas, allowing the for- 
mation of shocks within the settling gas itself; shocks 
can also be expected as a result of the interaction 
between infalling gas and the settling gas (Combes 
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et al. 1988: Athanassoula 1992). Interestingly, simi- 
lar results have been obtained in a different context, 
namely the disruption and orbital motion of stellar 
debris near a black hole (Kochanek 1994a; Lec ct al. 
1996). In such cases, the simulations clearly reveal 
collisions between material settling in the potential 
well and streams that arc still in the process of ac- 
cretion, with the resulting shocks contributing signif- 
icantly to a loss of angular momentum of the infalling 
streams. 

Preliminary three-dimensional hydrodynamic sim- 
ulations of interactions between a small gas-rich com- 
panion galaxy and a large elliptical (Koekemoer 1996), 
using parameters appropriate to those inferred for 
PKS 2356—61, indicate that infalling material is likely 
to interact directly with the settling gas for a compar- 
atively wide range of encounter geometries and ga- 
laxy/halo configurations, and that such interactions 
persist over timescales ~ 10'' — 10* yr. Direct con- 
nections between shock formation and tidal accre- 
tion have also been proposed elsewhere in the litera- 
ture, for example the peculiar spiral galaxy NGC 4438 
(Combes et al. 1998; Keel & Wehrle 1993; Kenney 
et al. 1995). Substantially more sophisticated mod- 
elling will be needed before more quantitative state- 
ments can be made concerning the properties of any 
shocks formed during the accretion process, in par- 
ticular detailed three-dimensional descriptions of the 
cooling and fragmentation processes occurring behind 
the shocks, as well as realistic models of the sizes 
and spatial distributions of gas clouds within the tidal 
streams. However, the results obtained to date sug- 
gest that shocks can play a significant role in galaxy 
interactions, with their energetic importance depend- 
ing primarily upon the velocities involved. 

Thus, the formation of shocks within the central 
EELR can be understood as being driven both by the 
dynamics of the settling gas itself, as well as by the 
additional infalling material. Since the timescale as- 
sociated with galaxy interactions is of the same order 
as the lifetime of the radio source, the EELR energet- 
ics can feasibly be powered in this scenario for much of 
the lifetime of the source. The energy supply for the 
EELR is ultimately related to the gravitational po- 
tential energy released by the infalling material, the 
kinetic energy associated with the initial relative ve- 
locity between the galaxies, and the total mass of gas 
accreted by the host galaxy. 



6. Conclusions 

We have presented results from detailed optical 
longslit spectroscopy of the EELR in the nearby pow- 
erful radio galaxy PKS 2356—61, providing complete 
three-dimensional spectroscopic datacubes that sam- 
ple the entire spatial distribution of line-emitting gas. 
We have carried out an extensive investigation of the 
observed physical properties of the gas, comparing 
these with results expected from the MAPPINGS II 
shock ionization code, and we have found that the 
excitation, kinematics and energetics of the gas are 
entirely consistent with a scenario involving auto- 
ionizing shocks as the dominant ionization mecha- 
nism. 

We find that the EELR is separated into two kine- 
matically and spatially distinct components: a cen- 
tral region, having a typical velocity width ^ 300 — 
400kms~^, and a tail region extending away from 
the centre, having much lower line-of-sight vc;locity 
widths (~ 120 — 150kms~^) and displaying more in- 
ternal velocity structure than the central region. We 
have also discovered a faint Ha-emitting galaxy lo- 
cated ~ 60kpc toward the north of the host galaxy 
of PKS 2356—61, exhibiting two elongated streams of 
emission apparently indicating a dynamical associa- 
tion between the two galaxies. 

The highly detailed nature of our data has allowed 
us to test in several cornplciinentary ways the feasi- 
bility of the proposition that the gas energetics may 
be powered by auto-ionizing shocks. Specifically, we 
have compared the EELR line-of-sight velocity widths 
with both the emission-line fluxes and excitation re- 
lationships expected from the MAPPINGS II mod- 
els, and we have also presented MAPPINGS II loci 
on diagrams of observed diagnostic line ratios plot- 
ted against one another. These tests repeatedly yield 
similar results: the properties of the central region 
are consistent with emission from shocks with intrin- 
sic velocities in the range ~ 400—550 kms~^, together 
with their associated precursors, while emission from 
the tail clouds are more indicative of "pure" precur- 
sor material, photoionized by shocks of these veloci- 
ties but not kinematically associated with them. This 
is further confirmed by an examination of the EELR 
energy budget, in which we find that a total covering 
fraction .Fcen ^ 0.6 for the precursor material in the 
central region is sufficient to account for the prop- 
erties of the central EELR, while at the same time 
providing agreement with the observed fluxes of the 
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photoionized clouds in the tail region. 

Although the alternative hypothesis of photoion- 
ization by radiation from an active nucleus cannot 
be completely ruled out in this object, we find that 
the radially uniform morphology of the central EELR 
cannot be readily reconciled with the anisotropic ion- 
ization "cone" anticipated in such a model. Further- 
more, nuclear photoionization requires additional ar- 
bitrarily free parameters, in particular the spectral 
shape and total flux of the ionizing continuum, to- 
gether with specific amounts of obscuration required 
to explain the fluxes in the tail clouds. On the other 
hand, the shock model can account for all the ob- 
served EELR properties by means of a single physical 
process, and allows direct observational verification 
of all the primary observable quantities — shock ve- 
locities, line fluxes and excitation — thereby greatly 
reducing the number of arbitrary parameters. 

The morphology of the EELR appears completely 
unrelated to the radio axis, and we find it unlikely 
that interactions between the radio plasma and the 
gas play a significant role in this source. Rather, 
we propose that the formation of the shocks is es- 
sentially related to gravitationally-driven kinetic en- 
ergy deposited by gas that is being accreted into the 
host galaxy from a companion object. While de- 
tailed three-dimensional radiative, non-equilibrium, 
magneto-hydrodynamic modelling would be required 
in order to determine the precise properties and life- 
times of shocks formed during such a process, the 
observed velocities correspond well to those required 
for gas accretion to take place during such inter- 
actions. Furthermore, the approximate interaction 
timescale in the case of PKS 2356—61 is in the range 
^ 10'' — 10^ yr, allowing for the possibility of a con- 
tinual supply of shocks from fresh infalling material 
during the lifetime of the radio source. 

Interactions between galaxies have long been con- 
sidered to play a role in triggering nuclear activity as 
a result of gas accretion, particularly in the context 
of powerful radio galaxies. Although the processes 
involved in such interactions are likely to be quite 
complex, the possibility that shocks can form as a re- 
sult, at least in some cases, and dominate the EELR 
energetics, is highly appealing in terms of accounting 
for all the observed properties of the extended line- 
emitting gas by a single physical process. Clearly, 
similarly detailed observations of other radio galaxies 
will be required in order to investigate the importance 
of this process in the general context of the formation 



and evolution of active galaxies. 
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Fig. 2. — (colour plate). Several views of the [OIIIJA5007 velocity datacube for PKS 2356—61, from dilFerent 
angles. The blue and green axes correspond to Right Ascension and Declination, with pixel dimensions of 1" and 
0'.'77 respectively. The pixel aspect ratio is such that the spatial image scale is the same along both axes. The 
red axis represents velocity, with increasing systemic velocity towards the right in all frames. The velocity scale is 
^ 45kms~^/pixel and the total velocity range is ^ 1800 kms^^. The colour scale corresponds to pixel flux (the 
brightest emission is coloured red, decreasing to yellow, green and blue). The flux scale is logarithmic, normalised 
to the rms noise level of the faintest detected emission (cto ~ 1 x 10~^^ ergs~^ cm~^ arcsec"^ A~^); values of Ictq, 
lOiTo: lOOiTo and ^ IOOOctq a-re indicated by dark blue, blue-green, yellow-orange and red, respectively. 



Fig. 4. — (colour plate). Spatial variation of emission-line flux ratios across PKS 2356—61: (a) [OIIIJA5007/H/5; 
(6) [NII]A6583/Ha; (c)_[SII]AA6716-|-6731/Ha; (d) [SII]AA6716/6731; (e) Ha/H/3. Pixel scales and image orien- 
tations are as for Fig. |l|. Colour indicates the ratio value, and ratios with relatively low S/N are shown with a 
correspondingly lower intensity. 
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